COOLING LINES AS PROBES OF THE FORMATION AND BUILDUP OF GALAXIES AND 

BLACK HOLES 

P.P. van der Werf 1 and M. Spaans 2 

1 Leiden Observatory, Leiden University, P.O. Box 9513, NL - 2300 RA Leiden, The Netherlands 
2 Kapteyn Astronomical Institute, Groningen University, P.O. Box 800, NL - 9700 AV, Groningen, The Netherlands 



Abstract 

We discuss the use of SPICA to study the cosmic his- 
tory of star formation and accretion by supermassive black 
holes. The cooling lines, in particular the high- J rotational 
lines of CO, provide a clear-cut and unique diagnostic for 
separating the contributions of star formation and AGN 
accretion to the total infrared luminosity of active, gas- 
rich galaxies. We briefly review existing efforts for study- 
ing high-J CO emission from galaxies at low and high 
redshift. We finally comment on the detectability of cool- 
ing radiation from primordial (very low metallicity) galax- 
ies containing an accreting supermassive black hole with 
SPICA/SAFARI. 
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galaxy. This is generally interpreted as evidence that the 
formation and growth of the SMBH is directly related to 
the formation process of the stellar population, e.g., in a 
violent burst of star formation. The picture of the simul- 
taneous build-up of a SMBH with an extreme burst of star 
formation provides a new context for the luminous and ul- 
traluminous infrared galaxies (LIRGs and ULIRGs), long 
suggested to be the birthplaces of SMBHs powering active 
galactic nuclei (AGNs). Studies of local ULIRGs demon- 
strate that they plausib l y evo lve into moderate mass field 
ellipticals (jGenzel et al.l . 120011 ). The giant ellipticals in the 
local universe would then be the result of similar but much 
more extreme events at higher rcdshifts. Indeed, the popu- 
lation of submillimetre galaxies (SMGs) provides the more 
luminous high redshift counterpart of the local LIRGs and 
ULIRGs and may be responsible for the formation of local 
massive spheroids and for g enerating QSO activity at high 



1. Introduction 

For several decades, our only direct probes of the high- 
redshift universe were QSOs and radio galaxies. While sig- 
nificant insight into the properties and evolution of these 
populations was obtained, almost nothing was known about 
the population of normal galaxies at high z. This has 
changed radically over the last decade. Deep imaging and 
spectroscopy from the ground and from space have re- 
vealed a rich and diverse population of high-z galaxies, 
and enormous progress has been made in establishing a 
complete inventory of the high-z universe, in terms of mass 
and metallicity budget, energy output and evolution over 
cosmic time. In addition, first investigations of the gene- 
sis of the fundamental scaling relations of galaxies (e.g., 
the fundamental plane and the Tully-Fisher relation) are 
being undertaken, and will provide insight into the origin 
of the galaxy population itself. 

One of the most remarkable scaling relations is the 
M m — a relation, i.e., the relation between the mass of a 
nuclear supermassive black hole (SMBH) and the velocity 
dispersion o f the spheroid (el li ptical or bulge) hosting the 
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Figure 1. CO emergent intensity distributions for a PDR 
and an XDR with identical impinging energy density, as 
described in Section [H Note that high-J CO emission is 
expect ed in both cases, but v astly more prominently in the 
XDR \Svaans & MeiierinA, \200A ). 



iFerrarese fc MerrittJ . 
l200ut h The fact that the black hole mass and total galaxy 
mass are closely related is truly remarkable, given that 
there is a factor of ~ 10 s between the AU-size Schwarzschild 
radius of the black hole and the kpc-size dimension of the 



2. Star formation history and accretion history 

Theoretically, a number of mechanisms h ave been pro- 
posed for establishing the M, — a relation (jMurrav et al.l . 
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20051: ICrotorJ . 2006), but attempts at observational tests 
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sions; (c) additional data for constraining local densities 
of these models are extremely cha llenging and have some- (density-tr acing molecules), column densities (isotopomeric 
times produced conflicting results ( Shields et al. . 20061 : AlexaiideB)atsid4 s of emission regions (high resolution ALMA 

imaging), as well as high-quality X-ray and infrared data. 
All of these will be in reach within the coming decade, 
and thus we may have some optimism that a reliable sep- 
aration of star formation and AGN as contributors to the 
total luminosity of galaxies (which has been elusive prob- 
lem for many years) can indeed be achieved with SPICA. 



SPICA provides a new window on this problem with 
its sensitivity in the far-infrared region. This region con- 
tains the key cooling transitions of photon-dominated re- 
gions (PDRs, heated by stellar UV radiation) and X-ray 
dominated regions (XDRs, heated by the X-ray radiation 
from an accreting SMBH). Advanced PDR/XDR models 
( Sp aans fc Meiierinkl . l2Q08h reveal the significantly differ- 
ent spectral patterns from XDRs and PDRs. In both cases 
cooling is provided by the principal fine-structure lines 
([C n] 158 jum, [O i] 63 and 145 /zm with in some regions of 
parameter space contributions from [Si II] and [Fe n] lines) 
as well as rotational lines of CO and rotational and vi- 
brational lines of H2. However, the flux distribution over 
the CO rotational ladder is vastly different for PDRs and 
XDRs. In the UV, a 10 eV photon penetrates a column 
of about 10 21 cm -2 (for solar metallicity) before being ab- 
sorbed by dust grains. UV photons also have the right en- 
ergy to dissociate CO between 11.2 en 13.6 eV. In contrast, 
a 1 keV X-ray photon traverses a column of 10 22 cm -2 
before absorption by hydrogen and metals occurs, since 
the absorption cross section scales as ~1 /energy 3 . Because 
of their high energy, neither X-rays nor the induced sec- 
ondary electrons destroy CO efficiently. Furthermore, UV 
photons cause photoelectric heating when electrons are 
liberated from (small) dust grains, with an efficiency of 
no more than 1%. X-rays on the other hand heat through 
ionization and (mostly) the subsequent Coulomb interac- 
tions of liberated very fast electrons, with an efficiency of 
about 10 — 50% (depending on the electron abundance). 
Consequently, PDRs, associated with star formation, con- 
tain CO layers that are relatively cool (10 — 100 K), while 
XDRs, associated with black hole accretion, are charac- 
terized by deeper layers of much warmer CO (20 — 500 K) . 

In Fig. Q] a comparison is presented between a PDR 
and an XDR for a density of n = 10 5 cm~ 3 and an il- 
luminating radiation field of 160 erg s -1 cm~ 2 , representa- 
tive of a region like the Orion Bar in the Milky Way (in 
the UV) or gas at a hundred parsecs from an accreting 
black hole powering a Seyfert nucleus with a luminosity 
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Figur e 2. Cooling budget in Mrk 231 iPavadovoulos et al. 

based on observations until CO J = 6—5 and a 
well- constrained model for the higher lines. The CO is 
modeled as a combination of a diffuse and a dense phase, 
with the dense phase containing most of the mass and 
emitting most of the cooling radiation, comparable in total 
luminosity to the [C II] line. 



(iMao et allfeood) and NGC 253 ([Bradford et all 120031 : iGusten et al 



20061 ). although for low-z galaxies they still do not probe 



3. State-of-the-art: low-z galaxies 

Reasonably detailed CO spectral line energy distributions 
of 10 44 erg s _1 (in X-rays). It should be noted that the are gradually be coming availab le, for instance for M82 
PDR case included here pushes the most extreme range 
expected in starburst galaxies, while X-ray fluxes as low 
as 5ergs _1 cm -2 , impinging on ~ 10 5 cm~ 3 gas are still 
capable of producing substantial CO emission at very high 
J CO levels. The CO ladder at high J levels thus consti- 
tutes a unique diagnostic of SMBH accretion activity in 
the SPIC A/SAFARI wavelength range. 

Real galaxies will of course produce a mixture of PDR 
and XDR emission. Separating these will require (a) a 
complete inventory of the cooling emission lines at high 



higher rotational levels than J = 7. The first high -quality 

CO flu x distribution for a ULIRG was presented bv lPapadopoulos et 
(|2007h for the nearby ULIRG Mrk 231. These authors con- 
clude that the CO ladder requires a two-component model: 
one component of fairly diffuse gas dominating the low- 
J lines and one component of dense gas dominating at 
J > 4. Remarkably, it was found that the dense compo- 



nent dominates the molecular gas mass as well as the total 
sensitivity; (b) accurate modeling with state-of-the-art PDR cooling radiation from Mrk 231. The total cooling by this 
and XDR codes, including radiative transfer in 3 dimen- component approaches that by the [C 11] line, implying a 
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totally different thermal balance than in for instance the 
Milky Way. The derived cooling budget is shown in Fig. [21 
This result also sheds new light on the discovery with the 
ISO satelite that in ULIRGs the [C ii]/Lir ratio is about a 
factor of 10 lower than in lower luminosity galaxies (where 
the [C n] 158 [im line typically contains 0.1 — 0.5% of Ljr). 
Physically, this situation can arise in dense PDRs (where 
the ionized carbon layer is thin and the molecular gas 
layer is close to the heating source), in XDRs, or in PDRs 
with an enhanced cosmic ray flux. These possibilities can 
be distinguished using the more comprehensive data that 
will be obtained by the Herschel Comprehensive ULIRG 
Emission Survey (HerCULES). This project will obtain 
full spectra of the CO ladder from J = 5 to 13 as well as 
the key cooling lines of [Cn] (158 /im) and [Oi] (63 and 
145 /im) for a flux-limited sample of 29 local LIRGs and 
ULIRGs. This project will be an important pathfinder for 
future studies with SPICA/SAFARI. 
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Figure 3. CO line luminosities in various high-redshift ob- 
jects compared to the Mrk231 template (bar diagram). 
Line luminosities are normalized to the Mrk231 values 
at eit her CO J=3— 2 or CO J=2— 1 t 'Pavadovoulos et all . 
20Qj) . 



4. State-of-the ART: HIGH-Z galaxies 

To date, detections of CO have been reported in about 50 
obj ects at z > 1, 38 of which are m entioned in the review 
bv ISolomon fc Vanden Bout! (|2005h . Of these 38, 22 are 
at redshifts z = 1.5 — 2.8, the time in cosmic evolution in 
which both the QSO and the SMG activity peak. Multi- 
line CO studies of the 12 brightest of these (8 QSOs and 
4 SMGs with no signs of AGN activity) are being carried 



out using the IRAM and APEX telescopes. These obser- 
vations typi cally cover CO rot ational lines from J = 3—2 
to J = 8-7 (jWeifi et al.l . l2005l ). A co mpilation, compared 



with t he Mrk231 flux distribution of iPapadopoulos et al 
(|2007h is shown in Fig. El 



The location of the peak of the CO flux distribution 
reveals different excitation conditions within the sample. 
For QSOs the peak typically occurs at the J = 6—5 or 7—6 
transition, indicative of a gas density of ~ 10 4 cm~ 3 and 
a kinetic temperature of ~ 50 K. For APM 0827+5255 the 
CO excitation is by higher by a large margin. Here the flux 
distribution peaks at the J = 10—9 transition. The SMGs 
in the sample tend to have a lower gas excitation than the 
QSOs with a peak of the CO ladder typically at the J = 
5—4 line. These CO ladders are similar to those observed 
in the central region of some local starburst galaxies (e.g., 
M82 and NGC 253, see references in Section[3]). The lowest 
CO excitation is observed toward the extremely red object 
HR 10 where the turnover of the flux distribution ahead ; 
occurs at the CO J = 4—3 line (jPapadopoulos fc Iviso 



2002). 



Although this set of CO lines constrains the proper- 
ties of the bulk of the CO emitting gas, surprising re- 
sults are obtained when the CO ladder is probed beyond 
the peak, as can be seen in the two QSOs at z ~ 4 
(APM 08279+5255 and PSS 2322+1944) where the obser- 
vations even include CO lines up to J = 11 — 10, in the rest- 
frame THz regime (jWeifi et all 120071) . The CO J = 9-8 
and 11-10 lines in the QSO host galaxy PSS 2322+1944 
(z = 4.1) show a clear flux excess compared to the single 
excitation model which describes the peak of the CO lad- 
der and require a second gas component with even higher 
CO excitation. This component may be related to a denser 
component (also indicated by HCN observations), or alter- 
natively to a hotter component heated by the AGN such 
as in APM 08278+5255, and can be modeled using the 
CO ladders shown in Fig. [TJ Interestingly, a similar CO 
ladder is observed towards the Galactic center, but there 
the high-excitation component is energetically much less 
important. 

In summary, observations of high-z sources allow ac- 
cess to diagnostic lines of CO which are very difficult to 
obtain from ground based observations for sources in the 
local universe. This leads to the somewhat paradoxical sit- 
uation, that the gas properties for some high-z redshift ob- 
jects are currently better constrained than those of most 
nearby infrared-luminous galaxies. The situation will be 
significantly improved by the HerCULES project, referred 
to in Section [3] Expanding such observations to lower lu- 
minosity galaxies, larger distances and higher J CO lines, 
in order to fully exploit the diagnostic use of the CO lad- 
der, will be a major task for SPICA/SAFARI. 
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Figure 4-. Rest frame spectral line distributions of high- J 
CO, H2 S(0) and S(l) (as well as lines from and 
HeH + ) are shown as functions of metallicity and X-ray 
flux for a density of 10 5 cm~ 3 . The impinging X-ray flux 
takes on values of 0.1, 1, 10 and 100 erg s^ 1 cmr" 2 ; and 



the metallicity value s of 1 , 
Svaans & Meiierink , 200 A ). 
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5. H 2 COOLING EMISSION FROM PRIMORDIAL GALAXIES 

A key question is what happens in galaxies of extremely 
low metallicity, where the [C 11] , [0 1] and CO cooling chan- 
nels are not available. In standard cosmological models, a 
significant population of dark matter halos collapses at 
redshifts larger than 5, aided by atomic (Hi and He 11 Ly- 
man a) and molecular (H2 and HD) cooling, and hosts th e 



2000). 



first stars and accreting black holes (jHaiman et al. 
The H2 rotational and ro-vibrational lines are accessi- 
ble to SAFARI and cool the (molecular) medium in sys- 
tems with metal l icities < 1 0~ 2 — 10~ 3 Solar. As s hown by 
Mizusawa et al.l (|2005l ) and lSpaans fc Meiierinkl (|2008l ), a 
number of these lines are expected to be detectable with 
JWST, ALMA and SPIC A/SAFARI, over a wide range of 
physical conditions. 

H2 cooling emission from population III star form- 
ing regions at z ~ 20 will be too faint for detection, 
but primordial star forming gas in z ~ 5 galaxies should 
be detectable with SAFARI due to the enhanced ioniza- 
tion and molecule forma tion caused by virialization shocks 
( Mizusawa et al. . 20051 ). Furthermore, the presence of an 
accreting black hole leads to the copious production of X- 
rays. In such X-ray irradiated gas, typical line intensities 
can be of the order of 0.1 — 1.0 erg s -1 cm~ 2 sr -1 for the 
pure rotational lines of H2 at 28, 17, etc. /im, which is 
as much as 10 — 30% of the black hole's acc retion-driven 
X-ray luminosity ( Spaans h Meiierinkl . 2008). Flux distri- 



butions for various parameters are shown in Fig. [4l which 
shows the dominant role of H2 at the lowest metallicities. 

Given the tight relation between supermassive black 
holes and the stellar spheroids that host them, it seems 
reasonable to use the extended Press-Schechter formalism 
and assume that (a) the M. — a holds at high redshift too, 
so the black hole mass is 2 x 10~ 4 of the total dark matter 
mass, that (b) black holes radiate at 10% of Eddington, 
yielding an X-ray dominated zone of 50 — 400 pc, and 
that (c) star formation and black hole growth are roughly 
simultaneous, with similar contributions to the luminos- 
ity. The latter implies star formation rates of 10 2 M Q yr _1 
within the central 2 kpc, quickly enriching the pristine gas, 
but the presence of a black hole alone is quite sufficient 
for SAFARI detections. 

For the ACDM concordance model, one then finds about 
10 of such sources per 2' x 2' field-of-view. The main con- 
tribution comes from z < 10 halos with masses > 10 9 M Q 
that yield H2 line fluxes of > 10~ 18 W m~ 2 , comfortably 
above the 2 x 10 -19 Wm~ 2 SAFARI detection threshold. 
Typically, at least two pure rotational H 2 lines (e.g., S(0), 
S(l), S(2), etc.) are detectable for a z ~ 5 source, allow- 
ing line confusion due to redshifted mid-infrared lines from 
far-infrared sources to be dealt with. SPICA/SAFARI tar- 
gets will likely be provided by JWST, which will discover 
many z > 5 galaxies. 
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